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Abstract 

In the SUSY GUT scenario, it is natural to assume the right-handed-neutrino 



> 

m 

p^J ' Major ana- mass scale to be 10 GeV. This will in principle lead, by the seesaw 

mechanism, to a Vt mass of order m^/(10^^ GeV) ~ 3 x 10~^ eV. This suggests 
\ that the solution of the solar-neutrino puzzle should be either the MSW effect in 

' ye~VT oscillations, with m^^ ~ 10~^eV^, or long-wavelength fe~f^ oscillations. 



with m^^ ~ 10"^'^ eV^. These solutions require unexpectedly large mixings of v^, 
with and f^, respectively. I suggest a variation of the Dimopoulos-Hall-Raby 
model for the fermion mass matrices which can accomodate such large mixings. 



Supersymmetric Grand Unified Theories (SUSY GUT's) have recently experienced 
^ ■ an upsurge in popularity, because they do not require the existence of an intermediate- 
■ energy breaking of the grand-unification group in order to obtain unification of the 
gauge couplings. In a SUSY GUT, the grand-unification gauge group is broken directly 
to the standard-model gauge group at Mq ~ 10^^ GeV, while supersymmetry is broken 
at a phenomenologically interesting scale Msusy ~ 10^ GeV [0. 

The seesaw mechanism for the suppression of the left-handed-neutrino masses 
can be incorporated in a SUSY GUT like SO(IO). That mechanism suggests that the 
mass of the Vr should be of order mf/Mj, where rrit is the top-quark mass and Mj 
is the scale of the breaking of the subgroup SU(4) [Q of SO(IO). This is because the 
right-handed-neutrino Majorana mass breaks SU(4), and should therefore be of order 
Mj] while SU(4) relates the neutrino Dirac mass matrix to the up-type-quark mass 
matrix, and therefore the z/^ Dirac mass should be of order rrit. The latter argument 
also suggests that the masses of the neutrinos should have a strong hierarchy. 

In the SUSY GUT scenario, Mj should be equal to Mq. We therefore guess that 
m^^ should turn out to be of order 3 x 10^'^ eV, and that m^^ <^ m^^ <^ m^,^. Such tiny 
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neutrino masses would be uninteresting from the points of view of direct laboratory 
measurements, of a solution to the dark-matter problem, and of a solution to the 
atmospheric-neutrino problem, but that should be accepted as a price of the SUSY 
GUT scenario ||]. 

The solar-neutrino problem can be solved by neutrino masses of these orders of 
magnitude in two different ways. One may use a Mikheyev-Smirnov-Wolfenstein 
(MSW) 0] resonant amplification of the Ve^^r oscillations in the sun. This requires 
10-^ eV^ and sin^(26'er) ~ 5 X 10 ^. One may alternatively use long- wavelength 



m 
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vacuum z/g-z/^ oscillations []^. This requires ml ~ 10 ^°eV^ and sin^(2^^e^) ~ 0.9 
i- 

Obtaining such large lepton-mixing angles is problematic in the context of SO(IO) 
models of the fermion mass matrices. In those models, the lepton mass matrices are 
related to the quark mass matrices. One would then expect the lepton-mixing matrix K 
to have matrix elements of the same order of magnitude as the corresponding matrix 
elements of the quark-mixing matrix V. This would lead to sin^(2^e^) ~ 10^^ and 
sin^(2^er) ~ 10~^, mixing angles much too small for solving the solar-neutrino puzzle. 

A nice model for the fermion mass matrices in a SUSY GUT is the Dimopoulos- 



Hall-Raby (DHR) model |10, 11|. This model has great predictive power in the charged- 
fermion sector: at the energy scale M/ = Mq-, it predicts that = m,-, ~ fri^/S, 

nid ~ Snie, \Vcb\ ~ ^Jmc/mt, and iKib/Kbl ~ ^mu/rric- The top-quark mass must be 
quite high, ~ 180 GeV, both in order to get \ Vcb\ small enough, and in order to obtain 
mft = rrir at the scale Mj, by using the fixed-point structure of the renormalization- 
group equations (RGE). 



The DHR model has been extended to the neutrino sector 12 and has encountered 



there the problem mentioned above. Indeed, in order to explain the solar-neutrino puz- 
zle, DHR have had to assume that it is the z/g-t'^ oscillations which are MSW-enhanced. 
That requires Mj ~ 10^^ GeV, two orders of magnitude smaller than Mq- This is incon- 
sistent with the basic philosophy of having a SUSY GUT without intermediate-energy 
symmetry breakings. 

The purpose of this Brief Report is to suggest a modification of the DHR scheme, 
which keeps some of its predictive power in the charged- fermion sector intact, while 
solving the solar-neutrino problem with Mj = Mq- I have explored various modifi- 
cations of the DHR scheme, and have found that they all predict too small a lepton 
mixing to be able to solve the solar-neutrino problem in a consistent fashion. The only 
exception that I have found is presented in this Brief Report. 

I suggest that at Mj = Mq the fermion mass matrices are 



Mfl = a 6 , (1) 

V d c / 

/ a -3d \ 

Me = \ a -3b , (2) 



2 



Mjj 



M, 



uN 



M, 



NN 



f 








s 


s 


n 


-3f 





n 


n 





S 





d 


b 





d 






X r . 



(3) 
(4) 
(5) 



These are, respectively, the mass matrix of the down-type quarks, of the charged lep- 
tons, of the up-type quarks, the Dirac mass matrix of the neutrinos, and the Majorana 
mass matrix of the right-handed neutrinos. The effective Majorana mass matrix of the 
left-handed neutrinos is given, in the seesaw approximation, by 
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a, b, c. 



d, /, s, n and r are complex numbers, a originates in the Yukawa couplings of 
the 10^ of SO(IO). c and n originate in the Yukawa couplings of the 10^, a different 
10. s originates in the Yukawa couplings of the 10^. b and d originate in the Yukawa 
couplings of the 126^. / originates in the Yukawa couplings of the 126^. The 126^ has 
a standard-model-breaking vacuum expectation value (VEV) contributing to the mass 
matrices of the down-type quarks and of the charged leptons, and a standard-model- 
invariant VEV leading to M^r^r. Similarly, the 10^ has VEV's in two standard-model- 
breaking directions, one VEV contributing to M^) and Me-, another one contributing 
to Mu and M^tv- 

The parameter r is the ratio between the standard-model-invariant VEV of the 
126^, and its standard-model-breaking VEV which leads to a contribution to Md and 
to Me- The latter VEV cannot be larger than the effective VEV included in Me, and 
Mei which is, because below Mq we have the minimal supersymmetric standard model 
with only two higgs doublets, v cos (3 = (175 GeV) cos /?. I now rely on the fact that the 
present model makes predictions for the charged-fermion sector which are very similar 
to the predictions of the DHR model, and I rely on the DHR analysis |T0|, to conclude 
that sin/5 should in the present model be larger than 0.9, a more likely value being 
much closer to 1 |^, just as in the DHR model. This means that f cos/3 < 75 GeV. 
Therefore, from the philosophy that SU(4) should be broken at Mq ~ 10^^ GeV, I find 
\r\ ~ 10^^ 

Notice that this model can be enforced by a simple Zn symmetry, n > 5, on the 
Yukawa couplings. With tu" = 1, the 10^ transforms as cu"^, the 10^ transforms as 
uj~^, the 10'^ transforms as uj~^, the 126"'^ transforms as ci;~^, and the 126^ transforms 
as Lu~'^. The three lepton generations are in representations 16^'^'^ of SO(IO); the 16^ 
transforms as uj^, the 16^ transforms as cj^, and the 16'^ transforms as lu^. I consider 
the existence of such a symmetry very important. In general, assuming the presence of 
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arbitrarily-located "texture zeros" in the mass matrices is inconsistent from the field- 
theory point of view. Often, one may even prove that there is certainly no symmetry 
capable of enforcing those zeros. 

I now work out the form of the quark- and lepton-mixing matrices in this model. 
I perform the analysis at the energy Mg, at which energy the mass matrices are as in 
Eqs. |l] to |. The purpose of the analysis is to show that one gets at that energy scale 
values for the charged-fermion masses and for the quark-mixing matrix very similar to 
the ones in the DHR model. If this is so, the machinery of the RGE then runs those 
values down to the weak scale in just the same way as it does in the DHR model, and 
finally the agreement of both models with experiment is just as good. One does not 



need to run the mass matrix themselves, as DHR have done |Ty, |T^; it is equivalent, 
but easier, to perform the bi-diagonalizations at the scale Mq and then to run the 
RGE for the masses and for the mixing-matrix parameters down to the weak scale; 



this procedure allows one to work with a smaller set of differential equations ||13 
Let the orthogonal matrix Ru be such that 



R 



/ I/I 

|s| I i?;7 = diag(m„, -mc,mt) . 
\ I'^l I'^l 



(7) 



Clearly, |/| = |?t,| = rrit — rric, and \s\ = y^m^m^. Therefore, 
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is a matrix which only mixes the second and third generations. The unitary matrix 
Ud diagonalizes M^Mj^, after some phases have been removed from that matrix: 



U 



D 



( \(^? + W 

\ah\ 
\ \cd\ 



\ab\ \cd\ \ 

|ap + |6p ladle*"^ | f//) = diag(m^, m^, m^) , 
\ad\e-'^ |cp + M|2 



(9) 



where = arg[(a^c)/(6(i^)]. Then, the quark mixing matrix V is given by 



V = R^diag{l,e'',l)UD, 



(10) 



where 9 = a.Tg[{bdn) / (acs)]. I have used the fact that Ru only mixes the second and 
third generations in order to eliminate one phase from V, by means of a rephasing of 
the up-quark field. 

Similarly, in the lepton sector, let the unitary matrix R^ be such that 
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-re*'^''*^^''/''^Miag(mi,m2,m3) . 
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mi, m2 and are real and positive, with mi < m2 < m^, and they are the hght- 
neutrino masses. The unitary matrix Ue diagonahzes MeM\^, after some phases have 
been removed from that matrix: 

/ |a|2 + 9|rf|2 -3|a6| -3|cd| \ 
Ul\ -3|a6| |ap + 9|6p -?>\ad\e^^ W/i? = diag(m^, mj, m^) . (12) 
V -3|cd| -?,\ad\e-'^ \c\'^ + 9\d\^ ) 

Then, the lepton mixing matrix K is given by 

K = i?^diag(e^«,e'^l)t/s, (13) 

where ^ = arg[((i^s^)/(6c/n)]. 

We see that this model involves seven real numbers (excluding |r|) and three phases 
in the determination of the fermion masses and of the mixing matrices, as opposed to 
the DHR model, which involves only six real numbers and one phase. 

If d vanished, Mr, and Me would be the same as in the DHR model. Ud would 
only involve a first-second-generation mixing. Because Ru only involves in the present 
model a second-third-generation mixing (see Eq. |^), one concludes that a vanishing d 
would lead to a vanishing Vub, which is excluded by experiment. One may indeed show 
that in this model \Vuh\ is proportional to \d\: 

« t(!!!i + cosV,), (14) 

ms + md mgmb \mb J 

iK.r ^ (15) 

mi 

Thus, the Cabibbo angle agrees with experiment even if c? = 0, but 

|c?|2 ^ m^lKfeT ~ m.mrf, (16) 

which means that \d\ and \a\ are of the same order of magnitude. From Eq. |l^, the 
result for the Cabibbo angle is not altered by the presence of such a small d, and 
iKisI ~ \Jmd/ms. Ud remains dominated by a first-second-generation mixing, and 

therefore \Vcb\ comes mostly from Ru- Just as in the DHR model, \Vcb\ ~ ^Jmjmt (see 
Eq. and therefore the top-quark mass must be high. The smallness of \d\ also makes 
for it that the relationships between the down-type- quark masses and the charged- 
lepton masses are essentially unaltered from what they were in the DHR model. The 
main difference between the two models is that iKife/Kbl was predicted to be small in 
the DHR model, while in the present model it is largely arbitrary, for it is proportional 
to the arbitrary parameter \d\ (see Eq. p!5|). 

Now consider the mass of the heaviest neutrino. As ^ \nf jd\ 3> |s//c?|, we 

find ma ^ (mcmt)/(ms|r|) ~ 3 x 10"'^ eV for |r| ~ 10^^. This is in the ball-park to 
explain the solar-neutrino deficit via MSW-enhanced v^-Vr oscillations. Notice however 
that this approximate formula is different from the one that had been guessed in the 
beginning of this work. 
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Another interesting feature of the model is that the mixing angle of the electron 
and tau neutrinos increases when \ Vub\ decreases, while the naive analogy between the 
lepton- and the quark-mixing matrices would have suggested the opposite. Indeed, \Vub\ 
is proportional to \d\. But, when \d\ decreases, the (1,3) and (2,3) matrix elements 
of M,yu increase, as seen in Eq. ||; this leads to larger mixings between Ue and z/,-, and 
between i/^ and Ur (but not between and u^). 

For the numerical work, I have used the approximations suggested by Naculich |jT3| 



The supersymmetry-breaking scale Msusy is fixed at 170 GeV, approximately equal 
to the top-quark mass. The one-loop RGE of the supersymmetric standard model are 
used to evolve the masses and mixing angles from Mg down to Msusy- The gauge 
couplings are analytic functions oit = \yi{M / Msusy) / {^^t^"^)'- 



585 - 5287rt 



9/ N 407r , 

glit) 



687 + 16807rt 



They unify for tc = 71/(1127r), at which scale they have the value g'^itc) = 357r/219. 
The unification scale is Mq = Msusy ^^pi^QTr'^tc) ■ For |r| I take, in each particular 
case, the value 

, , Mg 34 /7l7r^ 



\r\ is therefore a function of (3. Notice that |r| might be larger than the value in Eq. [1^, 
for the standard-model-breaking VEV of the 126^ is unknown and is only bounded to 
be smaller than f cos/5. |r| might therefore very well be one order of magnitude larger 
than the value that I take for it, which would correspondingly suppress all the light- 
neutrino masses. At Msusy, the values of the quark- and charged-lepton masses are: 
TTir = 1.749 GeV, = 103.4 MeV, = 500.9 keV, mb = 2.89 GeV, m, = 81MeV, 
rrid = 4.1 MeV, rric = 672 MeV, and = 2.4 MeV, as obtained by Naculich from the 
QCD and QED running of those masses from the energy scale at which they are known 
up to Msusy- The quark- mixing parameters do not run significantly at energy scales 
below Msusy- These are the values that I have tried to fit, taking especial care to 
fit correctly the charged-lepton masses and the Cabibbo angle. The top-quark mass is 
rrit = 1.04:3 mt{MsusY), due to the QCD correction. 

Let us denote the squared eigenvalues of the Yukawa- coupling matrix of the up-type 
quarks [Mjj /{v sin /?)] by Ua, and the squared eigenvalues of the Yukawa-coupling matri- 
ces of the down-type quarks and of the charged leptons [Mr, / {v cos /?) and Me / {v cos /?) , 
respectively] by Di and Ei, respectively. The indices a and i take the values 1, 2 and 
3. These squared eigenvalues run with energy according to 



dUa 

~dr 

dPj 
dt 



Ua{au + huUa + 2ca/ Di\Va. 

i 

Di{aD + 6a/A + 2CMY.Ua\Va^f) , (19) 
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Table 1: Examples of sets of input values. The phases tp, 9 and ^ are in degrees; all 
other values, except tan /3, are in GeV. 



dt 



EAaE + bMEi 



The iKtip are the squared moduli of the quark-mixing-matrix elements, which also run 
with energy, following Jl4| 



1 d\Vj^ 
Cm dt 



2 IK. 



-u^ - A + E + T.Dj\Vaj\ 
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+ 2 E E [2Re (K.v^ft-v;*,-!/;, 



w ■? 

Di - Dj 

( UgPj 



+ 



(20) 



Here 14, O 



2Re {v^^Vp,V:jV;,) = 1 - \VJ^ - - |V;,f - \Vp,\^ + |KiV>,f + K,Vpi\^ (21) 

for /3 7^ a and j ^ i. The parameters a^/, a^, a^;, hu and CAf are model-dependent. In 
the super symmetric standard model, cm = 1 and hu = 6, and 

26 32 
at/ = -6(72'-y^7| + 6(t/i + f/2 + f/3), 

14 32 

az? = -—gl-Q9l-—9l + QiD^ + D2 + D,) + 2iE, + E2 + E,), (22) 
= -^(7^-6(72' + 6(A + A + A) + 2(A + A + A)- 



I use as input the up-type-quark masses and the values of |a|, \b\, |c|, \d\, and the 
phases ip, 9 and ^ at the energy Mq- Some typical fits are given in table 1. 

The neutrino masses and the lepton-mixing matrix are calculated at Mq and do 
not run with energy. For each of the fits given in Table 1, the values of the neutrino 
masses and of the relevant lepton-mixing parameters are given in Table 2. For the 
calculation of the neutrino masses, I have assumed |r| to be given by Eq. |1^. Also, 
sin2(26'e^) = 4:\Kif.K2^Ki^K2e\i with similar definitions for sin2(2^e^) and for sin2(26'^^). 

One sees that it is possible to obtain sufficiently large to allow for the solar- 
neutrino problem to be solved by MSW-enhanced Ve^^r oscillations. Also, sin2(2^, 



is quite large, and it may be as large as to allow a solution of the solar-neutrino puzzle 
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Table 2: Neutrino-masses and mixing parameters yielded by the three examples in 
Table 1, respectively. 



by means of lye^i^n vacuum oscillations. The neutrino masses are in general adequate 
for those explanations of the solar-neutrino deficit. One sees that m2/mi ~ 3 and 
7713/777-2 ~ 300. 

In the present model, a solution of the solar-neutrino problem purely by Pe-J^ii 
vacuum oscillations is possible. The mixing angle ^er can be made sufficiently small 
to render the MSW effect between i/g and i^t irrelevant. This is what happens in the 
second example in Tables 1 and 2. On the other hand, a solution of the solar- neutrino 
problem solely by the MSW effect between Ue and 1/^-, though possible if one just 
takes into account the masses and the mixing angle of those two neutrinos (see the 
first example in Tables 1 and 2), should not be attempted. This is because the mass 
difference ttt.^ — trii is in general in the right range to suppress the '^B neutrinos by 
means of long- wavelength vacuum oscillations, and the mixing angle is always so 
large that the effect of those vacuum oscillations cannot be neglected. It is possible to 
find examples (with high tan/9) in which — ml is so small as not to suppress the 

neutrino signal meaningfully by means of the vacuum oscillations (the pp neutrino 
signal will still be suppressed); but in those examples, because 777.3/777.2 ~ 300 in this 
model, 77i| — 771^ will also be too small for the MSW effect to be able to explain the 
solar- neutrino deficit. 

In general, a mixed solution of the solar-neutrino problem should be considered. 
An example which requires such a mixed interpretation of the solar-neutrino deficit is 
the third one in Tables 1 and 2. There, ml — m\ a; 10-5 eV^ and sin^(2^eT) ~ IQ-^ 
are such that the MSW effect provides a suppression of the solar-neutrino signal in 
all present experiments (Homestake, Kamiokande and gallium) equal to about 50% of 
the suppression which is actually observed. The other half of the observed suppression 
is attributed to the fg-^'^ vacuum oscillations, with m\ — mf a 1.1 X 10-1° eV^ and 
sin^(2^e^) 0.68. In general, in the present model one should consider both the 
MSW effect in the Pe^i^r oscillations, and the long-wavelength Ve-i^n oscillations, when 
attempting to explain the solar- neutrino deficit, because both effects will suppress the 
solar-neutrino signal significantly. 

In conclusion, I have argued in this Brief Report that in the SUSY GUT scenario the 
neutrino masses should be extremely small, in such a way that the MSW explanation of 
the solar-neutrino deficit cannot be through Ug-Ui^ oscillations, but can quite likely occur 
through i^e-J^T oscillations. Vacuum i^e-J^n oscillations are also a possible explanation 
of the solar-neutrino deficit in this scenario. As an illustration, I have constructed a 
mass-matrix model which yields predictions for the charged-fermion sector similar to 
the ones of the DHR model, but which is able to fit the relatively large mixing angles 
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needed in these alternative views of the solar-neutrino suppression. 

I thank Lincoln Wolfenstein for first calling my attention to the problem of the over- 
all scale of the neutrino masses in SUSY GUT models in general, and in the DHR model 
in particular. I also thank him for discussions, and for reading the manuscript. This 
work was supported by the United States Department of Energy, under the contract 
DE-FG02-91ER-40682. 
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